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We evaluate static properties and semileptonic decays for the ground state of doubly heavy 
H, S', H* and S7, fl , £7* baryons. Working in the framework of a nonrelativistic quark model, we 
solve the three-body problem by means of a variational ansatz made possible by heavy quark spin 
symmetry constraints. To check the dependence of our results on the inter-quark interaction we 
use five different quark-quark potentials that include a confining term plus Coulomb and hyperfine 
terms coming from one-gluon exchange. Our results for static properties (masses, charge radii and 
magnetic moments ) are, with a few exceptions for the magnetic moments, in good agreement with 
a previous Faddeev calculation. Our much simpler wave functions are used to evaluate semileptonic 
decays of doubly heavy H, H'(J = 1/2) and S7, S7'( J = 1/2) baryons. Our results for the decay widths 
are in good agreement with calculations done within a relativistic quark model in the quark-diquark 
approximation. 

PACS numbers: 12.39. Jh,12.40.Yx,13.30.Ce,14.20.Lq,14.20.Mr 



I. INTRODUCTION 

The subject of doubly heavy baryons has been attracting attention for a long time. Magnetic moments of doubly 
charmed baryons were evaluated back in the 70's by Lichtenberg [l[ within a nonrelativistic approach. The infinite 
heavy quark mass limit was already used in the 90's to relate the spectrum of doubly heavy baryons to the one of 
mesons with a single heavy quark Q , or to analyze their semileptonic decay @ . A factor of two error in the hyperfine 
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TABLE I: Quantum numbers of doubly heavy baryons analyzed in this study. S, J P are strangeness and the spin parity of the 
baryon, I is the isospin, and S% is the spin parity of the heavy degrees of freedom. I denotes a light u or d quark . 
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splittings of Ref. have been recently noticed by the potential nonrelativistic QCD (pNRQCD) calculation of Ref. 0]. 

On the experimental side the SELEX Collaboration claimed evidence for the S+ baryon, in the h+K~ir + and 
pD + K~ decay modes, with a mass of M H+ = 3519 ± 1 MeV/c 2 @. Those results were challenged by a theoretical 
analysis [6( which claimed the observed events by the SELEX Collaboration could be explained without the involvement 
of doubly charmed baryons. Other experimental collaborations like FOCUS [7], BABAR Q and BELLE Q have 
found no evidence for doubly charmed baryons so far. At present the S+ has only a one star status and it is not 
listed in the particle summary table [Iol |. 

In hadrons with a heavy quark and working in the infinite heavy quark mass limit the dynamics of the light 
degrees of freedom becomes independent of the heavy quark flavor and spin. This is known as heavy quark symmetry 
(HQS) [HLQlGSEil. This symmetry was developed into an effective theory (HQET) [l5| that allowed a systematic, 
order by order, evaluation of corrections in inverse powers of the heavy quark masses. Unfortunately ordinary HQS 
can not be applied directly to hadrons containing two heavy quarks as the kinetic energy term needed in those systems 
to regulate infrared divergences breaks heavy flavor symmetry (l6j . For those hadrons the symmetry that survives is 
heavy quark spin symmetry (HQSS) fl7l |. which amounts to the decoupling of the heavy quark spins in the infinite 
heavy quark mass limit. In that limit one can consider the total spin of the two heavy quark subsystem (Sh) to be 
well defined. In this work we shall assume this is a good approximation for the actual heavy quark masses. This 
approximation, which is the only one related to the infinite heavy quark mass limit that we shall use, will certainly 
simplify the solution of the baryon three-quark problem. Recently the authors of Ref. [ij have developed and effective 
theory (pNRQCD) suitable to describe baryons with two and three heavy quarks. 

Solving the three-body problem is not an easy task and here we shall do it by means of a variational approach. 
The approach, with obvious changes, was already applied with good results in the study of baryons with one heavy 
quark [18]. This method, that leads to simple and manageable wave functions, is made possible by the simplifications 
introduced in the problem by the fact that we can consider Sh to be well defined. We shall consider several simple 
phenomenological quark-quark interactions [llJUliHH which free parameters have been adjusted in the meson sector 
and are thus free of three-body ambiguities. The use of different interactions will allow us to estimate part of the 
theoretical uncertainties affecting our calculation. Uncertainties related to the nonrelativistic baryon states that we 
use are difficult to estimate. We are aware of the limitations of a non-relativistic approach to describe light quark 
physics. However, for the kind of study performed in this work these are likely not as relevant as in other contexts. We 
study the semileptonic decays of double heavy baryons, in which the light quark is merely an spectator (heavy-to-heavy 
transitions). Indeed the lack of a proper relativistic treatment did prevent us to study semileptonic transitions of the 
type b —> u of a greater phenomenological interest. One should notice however that at least part of the relativistic 
effects not explicitly taken into account are included in an effective way in the parameters of the interaction which 
had been fitted to experimental data. We think this explains why the nonrelativistic quark model is so successful 
phcnomenologically even in the presence of light quarks. 

Our simple variational calculation reproduces the results for static properties obtained in Ref. 21| by solving more 
involved Faddeev type equations. Our method has the advantage that we provide explicit and manageable wave 
functions that can be used to evaluate further observables. Static properties like masses and magnetic moments of 
doubly heavy baryons have also been studied in other models. Masses have been calculated in the relativistic quark 
model assuming a light quark heavy diquark structure [HI], the potential approach and sum rules of QCD [23], the 
nonperturbative QCD approach [24j, the Bethe-Salpeter equation applied to the light quark heavy diquark |25| . the 
nonrelativistic quark model with harmonic oscillator potential [2(| or with the use of QCD derived potentials [27l . [28| , 
the relativistic quasi-potential quark model [29(, with the use of the Feynman-Hellman theorem and semi-empirical 
mass formulas within the framework of a nonrelativistic constituent quark model [30j . in effective field theories 0, 
l3~i| . or in lattice nonrelativistic QCD [32j]. There are also lattice QCD determinations [H, [H], [H|. Similarly, 
magnetic moments have been evaluated in a nonrelativistic approach |l|, in the relativistic three-quark model [361 ] . 
the relativistic quark model using different forms of the relativistic kinematics [37| , in the skyrmion model [38[ , in the 
Dirac equation formalism [3t| , or using the MIT bag model [4(| ■ 

We shall further use our manageable wave functions to study semileptonic decays of doubly heavy , J = 1/2, 
baryons. We shall evaluate form factors, decay widths and angular asymmetry parameters. Previous calculations of 
semileptonic decay widths have been done in different relativistic quark model approaches [4l|, [S^, EH , or with the 
use of HQET 0. 

The paper is organized as follows. In Sect. [H] we study the Hamiltonian of the system (Subsect. Ill A| and briefly 
introduce the different inter-quark interactions used in this work (Subsect. Ill Bl) . The variational wave functions are 
discussed in Sect. IIIII In Sect. HVI we present results for the static properties: masses (Subsect. IIV A[) . charge densities 
and radii (Subsect. IIV B[) . and magnetic moments (Subsect. IIV Cp . Semileptonic decays are analyzed in Sect. El After 
the presentation of general formulas, in Subsect. IV Al we relate the form factors to matrix elements and show how the 
latter ones are evaluated within our model. In Subsect. IV Bl we present our results for the form factors, differential 
and total semileptonic decay widths, and angular asymmetry parameters. The findings of this work are summarized 



FIG. 1: Definition of different coordinates used throughout this work. 



in Sect. IVIl The paper also includes three appendices: in appendix [X] we analyze the infinite heavy quark mass 
limit of our radial wave functions. In appendix [B] we relate the form factors for semileptonic decay to the two basic 
integrals in terms of which all of them can be obtained. Finally, in appendix [C] we give explicit expressions for those 
basic integrals. 

In Table U we summarize the quantum numbers of the doubly heavy baryons considered in this study 1 . 

II. THREE BODY PROBLEM 
A. Intrinsic Hamiltonian 

In the Laboratory (LAB) frame (see Fig.[T]), the Hamiltonian (H) of the three quark (hi, I12, q, where hi, /12 = c, b 
and q = l(u, d), s) system reads: 

/ V72. \ 

+ v hlh2 + V hiq + V h2q (1) 

where mh 1 , rrih 2 , m q are the quark masses and the quark-quark interaction terms Vjk depend on the quark spin-flavor 
quantum numbers and the quark coordinates (x% 1 ,x% 2 , x q for the hi, h%, q quarks respectively). To separate the 

Center of Mass (CM) free motion, we go to the light quark frame (R, fi, f^), 

ft m hl x hl + m h2 x h2 + m q x q 

K = 

m hl + m h2 + m q 

fl = X hl - X q 

r 2 = x h2 - x q (2) 

where R and fi , r*2 are the CM position in the LAB frame and the relative positions of the hi , /12 heavy quarks with 
respect to the light quark q. The Hamiltonian now reads 

V 2 - 

H = -^+H [at (3) 
2M 




1 Note that the definitions of H(, c and H! are interchanged in some references, with H^c having Sh = and H^ c having Sh = 1. The same 
applies to fi(, c and Q', . In tables we always quote the results corresponding to the convention we use. 
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H int = M+J2 H^ + V hlh2 {ri~r 2 ,spin) 



Vi • V2 



=1,2 

2 



H s P = + y hjq{?h spin)i j = 1,2 (4) 

where M = + mh 2 + = (l/ m, fcj + l/rriq) 1 and Vj = d/dp^ j = 1,2. The intrinsic Hamiltonian iJ lnt 

describes the dynamics of the baryon. Apart from the sum of the quark masses M, it consists of the sum of two single 
particle Hamiltonian Hj P , each of them describing the dynamics of a heavy- light quark system, plus the heavy-heavy 

interaction term, including the Hughes-Eckart term (Vi • V2)- We will use a variational approach to solve it. 



B. Quark— Quark Interactions 



We have examined five different interactions, one suggested by Bhaduri and collaborators [19( (BD) and four 
suggested by B. Silvestre-Brac and C. Semay [H El (AL1, AL2, API y AP2). All of them contain a confinement 
term, plus Coulomb and hyperfine terms coming from one-gluon exchange, and differ from each other in the form 
factors used for the hyperfine terms, the power of the confining term or the use of a form factor in the one gluon 
exchange Coulomb potential. All free parameters in the potentials had been adjusted to reproduce the light (it, p, K, 
K* , e tc.) and heavy-light {D, D* , B, B*, etc.) meson spectra 2 . All details on the above interactions can be found in 
Refs. Ill S3, El • ' 

These interactions were also used in Ref. [2ll ] to obtain, within a Faddeev calculation, the spectrum and static 
properties of heavy baryons. Our simpler variational method will not only give equally good results for the observables 
analyzed in [11, but it will also provide us with easy to handle wave functions that can be used to evaluate other 
observables. 



III. VARIATIONAL WAVE FUNCTIONS 

For the above mentioned interactions, we have that both the total spin and the internal orbital angular momentum 
given as 

S = (a hl + a hs + cr q ) /2 

L = h + fa, with lj — — i fj x Vj, j — 1,2 (5) 

commute with the intrinsic Hamiltonian and are thus well defined. We are interested in the ground state of baryons 
with total angular momentum J = 1/2, 3/2 so that we can assume the orbital angular momentum of the baryons to be 
L = 0. This implies that the spatial wave function can only depend on the relative distances r\, r-i and ri2 = \v\ — ri\- 
Furthermore when the heavy quark mass is infinity (mh — > oo), the total spin of the heavy degrees of freedom, 
Sheavy = (&hi + &h 2 ) /2, commutes with the intrinsic Hamiltonian, since the spin-spin terms in the potentials vanish 
in this limit. We can then assume the spin of the heavy degrees of freedom to be well defined. 

With these simplifications we have used the following intrinsic wave functions in our variational approach 3 

• H/nh 2 , Vl hl h 2 ~type baryons: 

\S hlh2 ,n hlh2 ;J=±,Mj) = ^^Ms.MsM^^-lMs^^W^Ms^ 

M Sh M Sq 

x*^ 2 (n,r 2 ,r 12 ) (6) 

where Mj is the third component of the baryon total angular momentum while \h1h2', Sh, Ms h ) and \q; ^Ms q ) 
represent spin states of the h^h^ subsystem and the light quark respectively, {jxjij^mxmzm) is a Clebsch-Gordan 



2 To get the quark-quark interaction starting from a quark-antiquark one the usual V?? = V?? /2 prescription, coming from a color 

dependence (A are the Gell-Mann matrices) of the whole potential, has been assumed. 

3 We omit the antisymmetric color wave function and the plane wave for the center of mass motion which are common to all cases. 
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coefficient. For hi = h 2 we need ^^'h (ri, r 2 , ri 2 ) = ^ifc (r 2 , ^12) to guarantee a complete symmetry of the 
wave function under the exchange of the two heavy quarks. 

• 5 iU2> K^-tyPe baryons: 

|5/U 3 > K lh2 ;J = \,Mj) = 0\\\Ms h M S9 Mj) \hxh»; lM Sh ) ® |g; l -M Sq ) 

M Sh M Sq 

x<Af(n,r2,r 12 ) (7) 

Similarly to the case before for h\ = h 2 we need 4^ ^ n r 2 , r 12 ) = (r 2 , r l5 r 12 ). 

• S'^^, W hlh2 -type baryons: 

\Z' hlh2 ,ti hxh2 ;J=\,Mj) = l/nftajOO)®!?;^) 

x^K'fri.ra.na) (8) 

In this case /ii ^ /i 2 and we do not need the orbital part to have a definite symmetry under the exchange of the 
two quarks. 

The spatial wave functions \&(ri, r 2 , 7*12) in the above expressions will be determined by the variational principle: 
5(B\H mt \B) = 0. For simplicity, we shall assume a Jastrow-type functional form 4 : 



*fcU 2 (ri,r 2 ,ri 2 ) = NF B {r 12 ) (t> hiq (n) cj> h2q (r 2 ) (9) 
where TV is a constant, which is determined from normalization 




with being the cosine of the angle between the vectors r\ and r 2 (ri2 = ( r\ + r 2 — 2rir 2 ^) 1 / 2 ). 

The functions fyhiq an< i <Phiq Wl U be taken as the S— wave ground states fjifj) of the single particle Hamiltonians 
Hj P of Eq. (0| modified at large distances. 

4>h ]q {r 3 ) = (l + ajrj)<pj(rj), j = 1,2 (11) 

The heavy-heavy correlation function F B will be given by a linear combination of gaussians 5 

4 

F B (r 12 ) = J2 a ^ (ri2+dj)2 > a i = 1 (12) 

3=1 

The value of one of the Oj parameters can be absorbed into the normalization constant N, so that we fix ai = 1. The 
rest of the variational parameters are determined by the variational condition <5(i?|-ff mt |£>) = 0. 

Although it is not self-evident from the functional form assumed, our variational wave functions are consistent with 
the infinite heavy quark mass limit, reducing in that limit to the product of the internal wave function for the heavy 
diquark times the wave function for the relative motion of the light quark with respect to a pointlike heavy diquark. 
All this is discussed in appendix |A"1 



4 A similar form lead to very good results in the case of baryons with a single heavy quark 1 1 SH . 

5 Note that F B should vanish at large distances because of the confinement potential. The confinement potential is also responsible for 
the non-vanishing values of the parameters ctj, j = 1,2 in Eq. illl . 
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10197 
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VAR 
FAD [21] 


10260 
10267 


10249 
10246 


10259 
10258 
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10224 


10274 
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a bb 


VAR 


10236 


10219 


10245 


10219 


10235 




VAR 


10297 


10287 


10301 


10269 


10302 




VAR 
FAD [21] 


6919 
6916 


6912 
6913 


6933 
6928 


6917 
6907 


6936 
6934 


£lbc 


VAR 
FAD [211 


6986 
7003 


6986 
6996 


6990 
6996 


6969 
6971 


7013 
7023 


—6c 


VAR 


6948 


6942 


6957 


6944 


6965 


Ml 


VAR 


7009 


7010 


7011 


6994 


7033 


—6c 


VAR 


6986 


6981 


7000 


6987 


6993 


ML 


VAR 


7046 


7047 


7055 


7037 


7057 



TABLE II: Doubly heavy S and Q baryqns masses in MeV. VAR stands for the results of our variational calculation. FAD 
stands for the results obtained in Ref. [2l| using the same interquark interactions but within a Faddeev approach. 



IV. STATIC PROPERTIES 



A. Masses 



The mass of the baryon is simply given by the expectation value of the intrinsic Hamiltonian. Our results (VAR) 
appear in Table |TT] where we compare them with the ones obtained in Ref. 2l| with the use of the same inter-quark 
interactions but within a Faddeev approach (FAD). For that purpose we have eliminated from the latter a small 
three-body force contribution of the type V123 = constant /m^m^rriq that was also included in the evaluation of 
Ref. [2l[ . We will show their full results in the following tables. Whenever comparison is possible we find an excellent 
agreement between the two calculations. In some cases the variational masses are even lower than the Faddeev ones. 
Besides we give predictions for states that were not considered in the study of Ref. [2l[ . 

In Tables iDTl and [TVl we compare our results with other theoretical calculations 6 . Our central values correspond 
to the results obtained with the AL1 potential, while the errors quoted take into account the variation when using 
different potentials. The same presentation is used for the results obtained in Ref. (2l| for which we now show their 
full values including the contribution of the three-body force. All calculations give similar results that vary within 
a few per cent. From the experimental point of view the SELEX Collaboration || has recently measured the value 
of Mw cc . This experimental value is 100 MeV smaller that our result. On account of what has been said in the 
introduction, one should take this experimental value with due caution. Note also that in Ref. Q the systematic error 
is not given. There are also different lattice determinations for baryons with two equal heavy quarks [33l l34l. |35| . Our 
results agree within errors with the lattice data for baryons with two c quarks, while they are roughly 100 MeV below 
lattice results for doubly 6-quark baryons. The best overall agreement with lattice data available so far is achieved in 
the calculation of Ref. [301 ] where they use the Feynman-Hcllmann theorem and semiempirical mass formulas in the 
framework of a nonrelativistic quark model but without the use of an explicit Hamiltonian. Within full dynamical 
calculations, ours and the relativistic calculations in Refs. [13, HH, have the best overall agreement with lattice 



Note in Ref. [3j| the S bc , 

well defined, being for S ba , Cl bc and 1 for E' bc , Q' bc - They are thus linear combinations of our states. The different spin functions are 
related by 

/ 1 \J=l/2 1 / 1 \ .7=1/2 /o / \ J=i/2 

(\ q c;0)®\b;±)) = - (j&c; 0} ® \q; ->J - ^ (j&c; 1> ® \q; ->J 

/ 1 \ -7=1/2 73/ 1 \ J=l/2 If I \ J=l/2 

(jgc; 1>®|&;-)J = (j&c; 0} ® \q; ->J - - (j&c; 1} ® [,; -> J (13) 

In order to extract their predictions for the 3(, c , S^ c and Q(, c , fl' bc baryons with total spin of the two heavy quarks well defined, we 
have assumed that the above relations, but with coefficients square, are also valid for the masses. Note this may be incorrect as we are 
neglecting a possible non negligible interference contribution. 



7 



data. 

This work [21] [22] [23] [24] [25] [26] ^7] [28] [29] [30] [31] [32] 

3612 +iY 3607+ 24 3620 3480 3690 3740 3646 3524 3478 3660 3660 ± 70 3610 3588 ± 72 

H* c 3706+ 23 3727 3610 3860 3733 3548 3610 3810 3740 ± 80 3680 

E bb 101971" 10194lig 10202 10090 10160 10300 10093 10230 10340 ± 100 

Z* bb 10236l? 7 10237 10130 10340 10133 10280 10370 ± 100 

E bc 6919±7 7 6915lg 7 6933 6820 6960 7010 6820 6950 6965 ± 90 f 6840 ± 236 

E' hc 6948t6 7 6963 6850 7070 6850 7000 7065 ± 90+ 

Z* bc 6986ls 4 6980 6900 7100 6900 7020 7060 ± 90 



This work Exp. [5] Latt. [33] Latt. [34] Latt. [35] 
3612+ iY 3519 ± 1 3605 ± 23 3549 ± 95 

3706+ 23 3 6 85 ± 23 3641 ± 97 

10197±^ 10314 ± 47 

10236+^ 7 10333 ± 55 

TABLE III: First panel: doubly heavy 3 masses in MeV as obtained in different models. Our central values, and the ones of 
Ref. [2l]], have been evaluated with the AL1 potential. Second panel: we compare our results with the experimental value for 
M~, cc measured by the SELEX Collaboration [5| (Note the cautions that appear on this experimental mass in the introduction) , 
and lattice results from Refs. [33l. l34l, |35(| . Entries with t should be taken with due caution (see footnote 6). 



-bh 



This work [21] [22][2^[2J[2J^[26][28][29J [30] [31] [32] 



Q. cc 3702+ 41 3710lf 3778 3590 3860 3760 3749 3590 3760 3740 ± 70 3710 3698 ± 65 



2+ 4i 37lFf 

Q* c 3783+ 22 3872 3690 3900 3826 3690 3890 3820 ± 80 3760 

n bb 10260t3 4 10267+. 4 3 10359 10180 10340 10340 10180 10320 10370 ± 100 

Q,Z b 10297+4 10389 10200 10380 10200 10360 10400 ± 100 

Ohc 6986+ 2 7 7003^3° 7088 6910 7130 7050 6910 7050 7045 ± 90 f 6954 ± 225 

Q! bc 7009t 24 " 7116 6930 7110 6930 7090 7105 ± 90 f 

Q.* bc 704619 1 7130 6990 7130 6990 7110 7120 ± 90 



This work Latt. [33] Latt. [34] Latt. [35] 

HZ 3702+ 41 3733 ± 9+3 8 3 663 ±97 

Q* cc 3783+ 22 3801 ± 9± 3 34 3734 ± 98 

Q bb 10260+3 4 10365 ±40;i2lJ 6 
nt b 10297^28 !0383 ± 39;gto 2 

TABLE IV: Same as Table [TTT] for doubly heavy Q baryons. 



There are also independent determinations of mass splittings in lattice QCD [H, H3, H^], nonrelativistic lattice 
QCD 32] and pNRQCD 0]. In Tablc|V] we compare those results to the ones obtained in the present calculation and 
in other models. Our central results evaluated with the AL1 potential are larger than the ones obtained in lattice 
QCD [3~3L [3^ . [H| and lattice nonrelativistic QCD [32]] . The agreement is better when we use the BD potential of 
Ref. [19] for which we always get the lowest results. Similar results are obtained by the relativistic calculation of 
Ref. [22|], whereas for the relativistic calculations in Refs. |23l [2a . [29j ] and the nonrelativistic one in Re f. [ 281) the 
agreement with lattice QCD and nonrelativistic lattice QCD data worsens. As for the calculation in Ref. [30( we do 
not quote their results due to the large theoretical errors involved. 

B. Charge densities and radii 

The baryon charge density at the point P (coordinate vector r in the CM frame, see Fig. [1]) is given by: 
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This work [22] [23J [25] [28] [29] [4] 



[32] Latt. [33] Latt. [34] Latt. [35] 
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TABLE V: Mass splittings in MeV for doubly heavy S and Q baryons. Our central values have been obtained with the AL1 
potential. Entries with a ' should be taken with due caution (see footnote 6). 
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FIG. 2: Charge form factor of the H cc and atb baryons evaluated with the AL1 20, 2] 
potentials. We show the two possible charge states. Similar results are obtained for H* 



| (solid line) and BD [ll| (dashed line) 
and E* bb . 



where eh 1 ,h 2 ,q are the quark charges in proton charge units e, and from Fig. [T]we have 7 yh 1 = y q + r\, ijh 2 = y, q + ^2 
and y q = — (m^ri + mh 7 f<i) / M. Since our L = wave functions only depend on scalars (r 1; and rri) the charge 
density is spherically symmetric (Pe(r) = pf(\f |)). 
The charge form factor is defined as usual 



(15) 



and it only depends on \q |. Its value at q = gives the baryon charge in units of the proton charge. 
The charge mean square radii are defined 



(r 2 )l 



/r+oo 
d 3 r r 2 pf{r) = 4tt J dr 



r*P* 



(r) 



(16) 



In Figs. [H [3] and [3] we show the charge form factors for S cc , Sbb, f2 cc , £lf,b and Sh c , f2b c baryons. In each case similar 
results are obtained for the star and prime excitations. We show the calculations with both the AL1 potential of 
Refs. [lOiHII and the BD potential of Ref. [ill]. The differences between the two calculations are minor in most cases. 

In Table [VI] we show the charge mean square radii. With the exceptions of the and Sl° c , we find good agreement 
with the results obtained in Ref. [2l| within a Faddeev calculation. The possible presence of a Sh — contribution in 



the wave functions of Ref. [2l| could be the possible explanation for this discrepancy. We also compare with the results 
obtained, for a few states, in Ref. (37j with the use a relativistic quark model in the instant form. The agreement is 
bad in this case. 



7 There exists the obvious restriction mf ll yf ll + m^y^ 




FIG. 4: Charge form factor of the H& c and fl bc baryons evaluated with the AL1 (solid line) and BD (dashed line) potentials. 
For S bc baryons we show the two possible charge states. Similar results are obtained for S bc , S' bc and fi£ c , fl' bc 



C. Magnetic moments 

The orbital part of the magnetic moment is defined in terms of the velocities v of the quarks, with respect to the 
position of the CM, and it reads 

M S = J d 3 ri d 3 r 2 (tffc lha (ri,r2,ri 2 ))* j^M^i x ™ hl v hl ) z 

+ Sr^ 2 X mh ^ h ^ z + 2^^« X m ^ z ] *fih 2 ( r i' r 2,ri 2 ) (17) 

with 8 m/ ll w^ 1 = — iVi, mh 2 Vh 2 = — and m q v q = i (V\ + V 2 ) • Since our orbital wave function has L = 0, the 
orbital magnetic moment vanishes. The magnetic moment of the baryon is then entirely given by the spin contribution. 

(B; J, Mj - J\ ^(a hl ) z + ^(a h2 ) z + ^(a q ) z \B; J, Mj = J) (18) 
Zm hl 2m h2 Zm q 

Those matrix elements are trivially evaluated with the results 



2 ( e hl _e^2 



i -" l '"- > "3 [2m hl + 2m h2 2 2m, 



Note that the classical kinetic energy has a term on Vf ll ■ Vf l2 and then the operator m bjl Vf ll is not proportional to — i V Skl , but it is 
rather given by m hl v hl = (M - m hl )/M ■ {-i Vy hl ) - m h2 /M ■ (~i V y)l2 ) = (-i Vi). Similarly m h2 v h2 = {-i XJ 2 )- 



10 



This work 



[211 



[37] 



This work 



' — cc 

77+ + 

1 — 'cc 


-0.030.™ 
0.298±g;°i| 


-0.038^ 

n Q1 r+0.035 


0.1 
1.3 


77* + 

' — 'cc 

H* ++ 


-0 42+ 007 

U.U^:Z_ 01 g 

0.341+°°12, 

— 0.042 






CO 
' — bb 

"bb 


0.221+°-^ 

w — 0.025 
n 1 qq+0.014 
u.ioo_ 016 


0.242+° 1 ° 3 i 

^ — 0.027 

-0 14S+ - 006 
u.iy:O_ 01g 




"bb 


-0.142+°-°" 

— 0.018 






77*0 
— bb 


fl OQQ+0.035 
u.zoo_ 031 






770 

"be 
77+ 
"be 


_n niS7+° 00 6 
U - UO '-0.013 

0.279l° ;™ 


-0.0721°;°°? 
0.306+°;°o?? 




C'O 
"be 
77'+ 
"be 


-0.0651°;°!° 

n ooq+0.036 
u -' oo -0 .025 






77*0 
"be 
77* + 
"be 


-0 065+ 0010 

U.UDO_ 01g 
U.OUO_ 039 







[21] 



[371 



file 


niQ+u.uui 

U.UlO_ 002 


0.009_ .oo3 0.2 




0.0091°,;°°,* 




«» 


-0.0861°;°°? 


-0 090+ 007 
u.uyu_ 002 


"bb~ 


-0 092+ 011 

U.U»Z_ ooi 




fib'e 


—0.016+0.003 


-0 n2 c i + ° °° 2 

U.UZO_ 003 


Q/0 

"be 


-0 1 9+ 003 

U.U1»_ 003 




o*° 

i2 bc 


-0 021 + 004 

U.UZ1_ 002 





TABLE VI: Charge mean square radii in fm 2 for doubly heavy S and f2 baryons. 
have been evaluated with the AL1 potential. 



Our central values, and the ones of Ref. 12 lH 
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(19) 



In Tabic fVTIl we give our numerical results. Our central values, as the ones obtained in Ref. [21| within a Faddeev 
approach, have been evaluated with the use of the AL1 potential. When compared to the values obtained in Ref. [2l| 
we find very good agreement with just a few exceptions (S° c , 2+^, H° c ). The discrepancy for the latter baryons may 
come from a possible non negligible Sh = contribution to their wave functions in the calculation of Ref. HH. In 
our case we have fixed S/j = 1 which we think is a very good approximation since in the limit of infinite heavy quark 
masses the spin of the heavy quark degrees of freedom is well defined. 

We also compare our results to the ones obtained in Refs.p], [H, H3, [H, [H, H3] using different approaches. The 
differences between different calculations are in some cases large. Being L = a good approximation the magnetic 
moments are essentially determined by the spin contribution of the quarks. With mj, 3> m u , rrid, m s , the contribution 
from the b quarks is negligible compared to the one of the light quark. This is also true to a lesser extent for the c 
quark. 



V. SEMILEPTONIC DECAY 

In this section we shall use the wave functions obtained with the variational method to study different doubly 
13(1/2+) — > B'(l/2 + ) baryon semileptonic decays involving ai»->c transition at the quark level. 
The differential decay width reads 



dT = 8\V cb \ 2 m B >G't 



d 3 p' 



d 3 k 



d 3 k' 



(2ii) 3 2E' B , (2tt) 3 2E Pi (2tt) 3 2E i 



7 (2n) 4 6 4 (p -p'-k-k 1 ) C a ^(k, k')H a p(p,p') 



(20) 



where \ V c b\ is the modulus of the corresponding Cabibbo-Kobayashi-Maskawa matrix element, is the mass of the 
final baryon, Gf = 1.16637(1) x 10 -11 MeV _2 [lO( is the Fermi decay constant, p, p', k and k' are the four-momenta 
of the initial baryon, final baryon, final anti-neutrino and final lepton respectively, and C and Ti are the lepton and 
hadron tensors. 

The lepton tensor is given as 



W(k,k') = k'»k a + k la k» - g^k ■ k' + ie^ aP k' a ko 



(21) 
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TABLE VII: Magnetic moments, in nuclear magnetons (|e|/2m p , with m p the proton mass), of doubly heavy E and fl baryons. 
Our central values, and the ones of Ref. [2lj . have been evaluated with the AL1 potential. 



where we use the convention e 0123 = 
The hadron tensor is given as 



-1, =(+,-,-, - 

r(o) 7 p(/-75)* 6 (o) 



B,rp) (B',r'p' 



* (0) 7ct (/- 75 )* (0) B,rp) (22) 



with \B, r p) (\B' , r' p'}) representing the initial (final) baryon with three-momentum pip') and spin third component 
r (r'). The baryon states are normalized such that (r p \ r' p') = (2ir) 3 (E(p)/m) 6 rr i S 3 (p — p'). The hadron matrix 
elements can be parametrized in terms of six form factors as 



B',r' p' 



* (0) 7m (/-75)* 6 (0) 



B,r p 



?i'(p'){l» (Fi(w) - 75G1H) + u„ (F 2 (w) - 75G2M) 

+v' ll (F 3 (w)- l5 G 3 {w))}u*(p) (23) 



where u B ' B are dimensionless Dirac spinors, normalized as uu = 1, and = p^/ms {v' — p'a/ m B l ) is the four velocity 
of the initial B (final B') baryon. The form factors are functions of the velocity transfer w — v ■ v' or equivalently of 
the four momentum transfer (q = p — p') square q 2 = m 2 B + m 2 B , — 2raBm B iw. In the decay w ranges from w = 1, 
corresponding to zero recoil of the final baryon, to a maximum value given by w = w max = {m 2 B + m B ,) / (2mBiTiB') 
which depends on the transition. 
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Neglecting lepton masses, we have for the differential decay rates from transversely (IV) and longitudinally (Tl) 
polarized W's (the total width is T = T L + T T ) HU 



dr T _ G%\v cb \ 



2 



dw 12n 
dT L _ G 2 F \V cb \ 2 



,3 



,^w 2 -lq 2 {{w - l)|fiM| 2 + (w + l)|G!(u;)| 2 } 

,Vw 2 -i{(w - i)\t v {w)\ 2 + (w + i)|^ A H| 2 } 



dw 247r 3 

F v ' A {w) = [(m B ±m B ')F? A (w) + (l±w)(m B ,F?> A (w) + mBFl> A (wj)], 
Fj = Fj(w), Ff = Gj(w), j = 1,2,3 



(24) 



One can also evaluate the polar angle distribution [45 

d 2 r 3 /dr T dr z 



dwdcosfl 8 \ dw 



2-p j jl + 2a' cos (9 + a" cos 2 61 j (25) 



where 9 is the angle between k' and p' measured in the off-shell W rest frame, and a' and a" are asymmetry 
parameters given by 

G 2 F \V cb \ 2 3 q 2 {w 2 -l)F 1 {w)G 1 {w) 
a = — 



6tt 3 a dr T /dw + 2dr L /dw 

dr T /du;-2dr L /du; 
dr T /du> + 2dr L /du; 



(26) 



These asymmetry parameters are functions of the velocity transfer w and on averaging over w the numerators and 
denominators are integrated separately and thus we have 

V ' 6tt 3 T T l + 2R L/T ' V ' 1 + 2R L/T ' L/T T T V ' 



A. Form factors 



To obtain the form factors we have to evaluate the matrix elements 



B',r' p' 



* C (0)7pa-75)* b (0) 



B,r p 



(28) 



which in our model are given by 
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(29) 



where the suffix "NR" denotes our nonrelativistic states and the factors d E/m take into account the different 
normalization. We shall work in the initial baryon rest frame so that p = 0, p' — —q, and take q in the positive 
z direction. Furthermore we shall use the spectator approximation. Having all this in mind we have in momentum 
space 
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x J d 3 qi d 3 q 2 ($?' h2 (qi 



m h2 +m q _ _ m h2 . B _ 

■9) 12 + =7?) I *6h 2 (9l:92j 



M> 



E b (qi) V 5 c (<fi - q) 



r+ Sl (<fi -«)7^(/-75)w Sl (ft) 



(30) 



where (<fi, 92) (^^/i 2 (9ii 9*2)) is the Fourier transform of the coordinate space wave function ^ bh2 (7*1, r2, fis) 
(f"i, ?*2, ^12 )) with gi, g*2 being the conjugate momenta to the space variables f\, r 2 . The factor of two comes 
from the fact that: i) for be— baryon decays, the charm quark resulting from the b — > c transition could be either the 
particle 1 or the particle 2 in the final cc baryon, while ii) for 66— baryon decays, there exist two equal contributions 
resulting for the decay of each of the two bottom quarks of the initial baryon. 
The actual calculation is done in coordinate space where we have 
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(31) 



where I = —i Vi represents an internal momentum which is much smaller than the heavy quark masses mb, m c . On 
the other hand \q \ can be large 9 . Thus, to evaluate the above expression we have made use of an expansion in I, 
introduced in Ref. [46], where second order terms in I are neglected, while all orders in \q\ are kept. For instance 
E C {1 — q) is approximated by E c (l — q ) « E c (q) X (1 — I ■ q /E 2 (q)) with E c (q) = yjm 2 c + q 2 . 
The three vector and three axial form factors can be extracted from the set of equations 10 
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for the vector form factors and 
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(32) 
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9 At q 2 = one has \q | = (m 2 B — m 2 B , )/2ms which is kj ma/3 for the transitions under study. 
10 Remember q is in the 2 direction. Notice also that, for p = 0, w is just a function of \q\. 
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fl',1/2 -q * C (0) 7o 75* b (0) 5,1/2 0) = 



' E B '{-q) + m B > 
2m B > 



E B >(-q) + m B > 



-G l {\q\) + G 2 (\q\) 



E B >{-q] 
m B , 



G,{\q\) 



(33) 



for the axial ones. All the left hand side terms can be evaluated using Eq. (|3ip with the approximation mentioned 
above. 

For each transition there are only two different coordinate space integrals from which all different matrix elements 
can be evaluated. Those integrals are 



i B ' B m 



K B ' B {\q\ 



j3 j3 t — f q-r-2 —i 1- q-r\ 

an a r 2 e M ' e m> 
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^ I d 3 ri d 3 r 2 e l T^ q ~ 2 e 



% h-2 + m <j - - 
— q-ri 



tffjLfa.ra.ria) 1 • 3 ^f/Jn^m) (34) 



In appendixlBlwe relate the form factors to the integrals I B B {\q \) and /C s B (\q \) for the different Sh, S' h combinations, 
while in appendix [C] we give the actual expressions we use to evaluate those integrals. 



1. Current conservation 



In the limit m& = m c and for B' = B (and thus S/j = S' h ) vector current conservation provides a relation among 
the vector F 2 and F3 form factors, namely 

F 2 (w) = F 3 (w) (35) 

On the other hand, if the two quarks in the current had the same flavour, vector current conservation would fix the 
normalization of the zeroth component of the vector current matrix element at w = 1, since this just counts the 
number of heavy quarks, two, in that case. For quarks with distinct flavours, but still in the limit m^, = m c , the 
forward vector matrix element has an extra Clebsch-Gordan factor of l/y/2, and thus we have 

F 1 (l) + F 2 (l)+F 3 (l) = V2 (36) 
In this limiting situation the integrals I BB (\q |) and IC BB (\q |) are related by 11 

K B Bm = m^^ xBBm (37) 

Besides one has that 1 BB (0) = 1. 

Using now the relations in Eq. (|B4[) in appendix 151 we see that our model satisfies the constraint in Eq. (f3"6"]) exactly. 
On the other hand we have a small violation of vector current conservation due to binding effects. For instance, and 
again using the relations in Eq. (|B4[) , we obtain for w = 1 

F 2 (1)=F 3 (1) + V2(1-^) (38) 

which shows that current conservation is violated by a term proportional to the binding energy of the baryon divided 
by the sum of the masses of its constituents. This violation disappears in the infinite heavy quark mass limit. This 
deficiency is avoided in most calculations by the neglect of binding effects between the heavy diquark and the light 
quark 12 . Improvements on vector current conservation would require at minimum the introduction of two-body 
currents [47[, going thus beyond the spectator approximation, that we have not considered in this analysis. 



One just has to integrate by parts in the KL BB {\q\) expression 
In Refs. [4ll . |4;J the currents are constructed at the diquark level. At the baryon level, vector current is conserved because light 
quark-heavy diquark binding effects are neglected. The calculation in Ref. [44| | uses the infinite heavy quark mass limit thus canceling 
binding effects. 
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FIG. 5: Vector F±, F2, F3 and axial G%, G2, G3 form factors for doubly H(J = 1/2) baryons decays evaluated with the AL1 
potential. The results for £l(J = 1/2) baryons (not shown) are very much the same. 



Note also that, for transitions that do not conserve the spin of the heavy quark subsystem Sh (i.e. S^, — > z^ c lvi) 
we have in the mj, = m c limit and at zero recoil that 

Fi(l)+F 2 (l)+F 3 (l) =0 (39) 

due to the orthogonality of the initial and final baryon wave-functions. 



B. Results 



In Fig. [5] we show the form factors for 5 decays evaluated with the AL1 potential. Variations when using a different 
potential are at the level a few per cent at most. The results for doubly heavy fl decays are almost identical to the 
corresponding ones for doubly heavy 3 decays. The fact that we have two heavy quarks and that the light one acts 
as a spectator makes the results almost independent of the light quark mass. 

In Fig. [6] we show now our results for the differential cffV/dw, dT^/div and dT/dw decay widths evaluated with 
the AL1 and BD potentials. The differences between the results obtained with the two inter-quark interactions could 
reach 30% for some transitions and for some regions of w. As a consequence of the apparent SU(3) symmetry in the 
form factors we also find that the results for doubly heavy S and f2 decays (not shown) are very close to each other. 
This apparent SU(3) symmetry goes over to the integrated decay widths and asymmetry parameters. 

In Table IVTlTl we give our results for the semileptonic decay width (transverse Ft, longitudinal Tl and total T) 
for the different processes under study. Our central values have been evaluated using the AL1 potential while the 
errors show the variations when changing the interaction. The biggest variations appear for the BD potential for 
which one obtains results which are larger by 7 ~ 12%. In Table ITXl we compare our results to the ones calculated 
in different models. For that purpose we need a value for \V c b\ for which we take \V c b\ = 0.0413. Our results are in 
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FIG. 6: dF/dw, dTh/dw and dTr/dw semileptonic decay widths in units of |V C (,| 2 ■ 10 -11 GeV, for doubly S(J = 1/2) baryons 
decays. Solid line, long-dashed line and short-dashed line: dT/dw, dTh/dw and dTr/dw evaluated with the ALI potential; 
dotted line, long-dashed dotted line and short-dashed dotted line: dT/dw, dTi,/dw and dTr/dw evaluated with the BD 
potential. The results for Q decay (not shown) are very similar. 



reasonable agreement with the ones in Ref. [41( where they use a relativistic quark model evaluated in the quark- 
diquark approximation 13 . The result for T(Sft c — > E cc ) obtained in Ref. [44[ using HQET is a factor of 4 larger than 
ours 14 . We can not directly compare our result for the latter transition with the one published in Ref. [13], the reason 
being that in this reference the S;, c state is defined differently, with the c and the light quark coupled to total spin 1. 
Had they defined the state with the be pair coupled to spin 1 they would had obtained roughly a factor of 4 larger 
result 50], in reasonable agreement with our calculation or the one by Ebert et al. 15 . Finally in Ref. [43[, where they 
use the Bethe-Salpeter equation applied to a quark-diquark system, they obtain much larger results for all transitions. 

In Table IXl we compile our results for the average angular asymmetries a' and cv", as well as the Rl/t ratio, 
introduced in Eq. pT)) . The central values have been obtained with the ALI potential. Being all quantities ratios the 
variation when changing the inter-quark interaction are in most cases small. 



13 Note we have divided by a factor 2 the results originally published in Ref. |4ll as there, and as we initially did, the authors forgot a 
symmetry factor 1 / y/2 for the case of diquarks with two equal quarks [48| . 

14 Note we have multiplied by a factor 2 the result originally published in Ref. [3 as the author overlooked a normalization factor y/2 l4tl . 

15 Note the actual physical states would be an admixture of what we call H;, c H^ c , and similarly in the f2 case. A measurement of the decay 
width could give information on the admixtures. 
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TABLE VIII: Semileptonic decay widths in units of |V c t| 2 ■ 1CP 11 GeV. Ft and Fl stand for the transverse and longitudinal 
contributions to the width F. The central values have been obtained with the AL1 potential. I stands for a light charged lepton, 
I = e, n 



This work [41] f [42]* [431 [_41 § This work [41] * 



F{E bb 


— > SbcWi) 


1 Q2+ - 25 

l- az -0.05 


1.63 


28.5 


F(n bb 


— > Qbc Ivi) 


2 14+ ' 20 

z - x ^-0.02 


1.70 


F{E bc 


— > Hcc Wi) 


+0.26 
z - d ' -0.03 


2.30 


0.79 8.93 8.0 


F(n bc 


— > n cc 


2,59+0.20 


2.48 


F{E bb 


-> Kcln) 


i 06+ 013 

l- uo -0.03 


0.82 


4.28 


T(Sl bb 




1.16+ 013 


0.83 


r(E' bc 


— > H cc Ivi) 


1 S6+ 010 


0.88 


7.76 


mi 




x 3 6 +o.o9 


0.95 



TABLE IX: Semileptonic decay widths in units of 10 14 GeV. We have used a value \V c b\ = 0.0413. I stands for a light charged 
lepton, I = e, fi. For results with t, J and § see text for details. 

VI. SUMMARY 

We have evaluated static properties and semileptonic decays for the ground state of doubly heavy 2 and fl baryons. 
The calculations have been done in the framework of a nonrelativistic quark model with the use of five different inter- 
quark interactions. The use of different quark-quark potentials allows us to obtain an estimation of the theoretical 
uncertainties related to the quark-quark interaction. In order to build our wave functions we have made use of the 
constraints imposed by the infinite heavy quark mass limit. In this limit the spin-spin interactions vanish and the total 
spin of the two heavy quarks is well defined. With this approximation we have used a simple variational approach, 
with Jastrow type orbital wave functions, to solve the involved three-body problem. 

Among the static properties, our results for the masses are in very good agreement with previous results obtained 
with the same inter-quark interactions but within a more complicated Faddeev approach In some cases we 

even get lower, and thus better, masses 16 . We have calculated mass densities and charge densities (charge form 
factors) finding that the corresponding mean square radii are again in good agreement with the Faddeev calculation 
of Ref. [2l[. We have also evaluated magnetic moments. Being the total orbital angular momentum of the baryon 
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TABLE X: Averaged values of the asymmetry parameters a and a" evaluated as indicated in Eq. (|27p . We also show the ratio 
Rl/t = Tl/Ft- The central values have been obtained with the AL1 potential. I stands for a light charged lepton, I — e, [i 



For a given Hamiltonian a variational mass is an upper bound of the true ground state mass. 
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L = 0, the magnetic moments come from the spin contributions alone. With the exception of S° c , c^, and £1® C we agree 
perfectly with the Faddeev calculation in Ref. [2l[. For the magnetic moments of 5° c , and f2° c the discrepancies 
between the two calculation are very large. The origin might be attributed to the presence of a non-negligible Sh = 
component in the wave functions of Ref. plj . In our case we have Sy L — 1 which we think is a good approximation 
based on the infinite heavy quark mass limit. This assertion seems to be corroborated by the results obtained in the 
relativistic calculation of Ref. (36[, at least for the S° c and f2° c cases. 

We have used our wave functions to study the scmilcptonic decay of doubly 5(J = 1/2) and S1(J = 1/2) baryons. 
We have worked in the spectator approximation with one-body currents alone. In the = m c case and for B = B' 
baryons we have checked that our model satisfies baryon number conservation. On the other hand we have a small 
vector current violation by an amount given by the binding energy over the mass of the baryon, violation that 
disappears in the infinite heavy quark mass limit. Small vector current violations due to binding effects seems to 
be present in most calculations to date. With this model we have evaluated form factors, asymmetry parameters, 
differential decay widths and total decay widths. Our results for the latter are in reasonable agreement with the ones 
obtained in Ref. (4ll | using a relativistic quark model in the quark-diquark approximation, while they are much smaller 
than the ones obtained in Ref. [43| by means of the Bethe-Salpeter equation applied to a quark-diquark system. 

For the weak form factors the results exhibit an apparent SU(3) symmetry when going from 5 to Q baryons. This 
is due to the fact that we have two heavy quarks and the light one acts as a spectator in the weak transition. This 
apparent symmetry appears also in the decay widths and asymmetry parameters. On the other hand SU(3) violating 
effects are clearly visible in some static quantities like the charge form factors and radii, and the magnetic moments, 
that depend strongly on the light quark charge and/or mass. 
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APPENDIX A: INFINITE HEAVY QUARK MASS LIMIT OF OUR VARIATIONAL WAVE FUNCTIONS 

In the infinite heavy quark mass limit the wave function of the system should look like the one for a "meson" 
composed of a light quark and a heavy diquark. The two heavy quarks bind into a 3 color source diquark that to the 
light degrees of freedom appears to be pointlike [3j. In our model the pointlike nature of the heavy diquark comes 
about through the one-gluon exchange Coulomb potential which binds the two heavy quarks into a distance 1 

1 m hl m h2 

rh t h 2 oc ; fj, hlh2 = ■ (Al) 

M/ufca m h! + m h2 

that tends to zero if both quark masses go to infinity. 

For our wave functions we can define the probability distribution Phim for the two heavy quarks to be found at a 
distance rh 1 h 2 

Phlh2 



(r hlh2 ) = j d 6 n J d 6 r 2 5(r 12 -r hlh2 ) \^ lh2 (r 2 ,r 2 ,r 12 )\ (A2) 

In Fig. [7J we show the Phim probability distributions for the S cc , S;, c , and £l cc , Of, c , &bb baryons evaluated for 
the AL1 potential. We see how the maximum moves, as expected, to lower rh 1 h 2 values as the quark masses increase. 

On the other hand the relative wave function for the light quark with respect to the heavy quark subsystem should 
tend to the one of a light quark relative to a pointlike diquark. This limit is not evident in the coordinates we work 
as one has first to separate the diquark internal orbital wave function from the total one. To see that this limiting 



17 The relation can only be approximate due to confinement and the interaction with the light quark. 
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FIG. 7: Ph 1 h 2 probabilities (see Eq. (lA2[) ) for the H cc , Ej, c , Sm> and £} cc , £}f, c , f7&b baryons evaluated with the AL1 potential. 



situation is in fact reached within our variational ansatz let us do the following: introduce the new set of coordinates 

m hl x hl + m h2 x h2 + m q Xq 



R 

r„ 



min + m h2 + m q 

•"hi *^/i2 

m hl x hl + m h . 2 x h2 



m hl + m h2 



m hl ri + m h2 r 2 
m hl + m h2 



in terms of which the Hamiltonian now reads 



where 



V 2 5 



H = -^ + H int 
2M 

H mt = M + H hlh2 +H qhl 



qh\}i2 



(A3) 



(A4) 



Hhlh 2 
Hqhih 2 



Vis 



if 

2 \m hl + m h 



+ Vfuhaifti, spin) 



1 



V 9 +V hiq (f q + 
+Vh 2q {r q - 



m h2 



m hl + m/,; 
mh t + m h 



■r*i2, spin) 
-ri 2) spin) 



(A5) 



with Vi2= 9/dff 12 and V g = d/dp q . Defining now 



£1 qh 1 h 2 



1 



1 



2 V 77 ^! + m fc 



one would have 



int 



h\}i2 



+ 



H qhih 2 



V„ +Vft ig (f g , spin) + V/ i2g (f,j, spin) 



qh\h 2 



n qh x h 2 



(A6) 



(A7) 



Hfnh 2 is the Hamiltonian for the relative movement of the two heavy quarks while H® hih is the Hamiltonian for the 
relative movement of the light quark with respect to a pointlike heavy diquark where the two heavy quarks are located 
in their center of mass. Both movements are coupled through the term (Jl q hih 2 — H^ hih ^). If the heavy quark masses 
get arbitrarily large the average ri2 value tends to zero and thus one can neglect the effect of Hqh^ — Hq\ x h 2 ■ ^ n that 
limiting situation the light and heavy quark degrees of freedom decouple completely and the internal Hamiltonian 
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reduces, as it should, to the sum of the Hamiltonians Ht ll h 2 i corresponding to the relative movement of the two heavy 
quarks, and H® hifl2 , corresponding to the relative movement of the light quark with respect to the pointlike heavy 
diquark subsystem. 

As to the wave function it should reduce in that limit to the product ^h 1 h- 2 { r i2) • ^qh 1 h 2 i r q)^ being 
<&h 1 h 2 ( r i2), <&qh 1 h 2 ( r <i') tne g roun( i state wave functions for Hh^, H® hlfl2 respectively. To check that we reach that 
limiting situation we have evaluated the projection V of our variational wave functions onto $hihs ( r i2) ■ ® q hih 2 ( r ?) 
evaluated with the actual heavy quark masses. This projection is given by 18 

V=J d 3 n J d 3 r 2 (^ lh2 (r u r 2 ,r 12 y $ hlh2 (r 12 ) $ qhlfl2 (r q ) ( A8 ) 

and the values for \V\ 2 that we obtain for the S cc , Sf, c , Ey, and fl cc , flbc, ^bb baryons using the AL1 potential are 
given in Table IXll We see how \V\ 2 increases with increasing heavy quark masses indicating that for very high heavy 
quark masses the total wave function tends to the one of a light quark relative to a pointlike diquark times the diquark 
internal wave function. In summary, our variational wave functions respect the infinite heavy quark mass limit. 



^cc ^bc ^bb 






\Vf 0.974 0.975 0.991 


0.959 


0.966 0.984 



TABLE XI: Absolute value square of the V projection coefficient defined in Eq. (|A8f) 



We note by passing that the product ^hih 2 { r v^) ' ^ q hih 2 ( r <?) corrected by a correlation function in the variable 
T\ 2 ■ f q would have been a good variational orbital wave function. For instance, a calculation using the AL1 potential 
and an orbital wave function given just by that product $hih 2 (712) • ® q h 1 h 2 ( r «) §i ves f° r the expectation value of H lnt 
for the S cc , E bc , S 66 baryons 

(H int )\t L1 = 3640 MeV; (iP nt )|^ L1 = 6943 MeV; (H in *)|£ il = 10198 MeV (A9) 

which are respectively 28 MeV, 24 MeV, and 1 MeV larger, and therefore worse (see footnote 13), than our best values 
in Table [ill The correlation function would clearly be needed, being its role more important for the "cc" and "6c" 
systems. The improvement in going from a "cc" to a "6c system" is not as good as one would naively expect, the 
reason being that, considering f\ 2 to be a small quantity, Hqh^ — H® hih2 is first order in r\ 2 for the "6c" system 
while, for symmetry reasons, it is necessarily of second order for the "cc" and "66" ones. In any case one sees how 
the expectation values obtained with this simple ansatz get closer to our best values in Table [ill as the heavy quark 
masses increase. 



APPENDIX B: FORM FACTORS IN TERMS OF THE X B B ( \q\) AND K. B B {\q\) INTEGRALS 

In this appendix we relate the vector F\, F 2l F3 and axial G\, G 2 , G3 form factors, that we evaluate in the center 
of mass of the decaying baryon, to the integrals T B B (\q |), K. B B (\q |) defined in Eq. (|3^)l . To simplify the expressions 
it is convenient to introduce 



S/ Un E B ,(-q)+m B > / 2E c (q) 



n E B ,{-q)+m B , / 2E c (q) . 



• Cases S h = 1, S' h = or S h = 0, S' h = 1: 

Fl(\q\) _ [2 ( l B ' B {\q\) IC B ' B (\q\) ( m c 



(Bl) 



E B ,(-q)+m B , \3\E c (q) + m c 2 \E 2 (q) m b 



Note r q can be expressed in terms of n t ra,ri3 and that dP , r\d 3 r2 = d 3 ri2d 3 r q 
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Fl(\q\) | F 3 (|g|) = 

E B ,(-q)+m B > m B ' 

h{\q\) + F 2 (\q\)+ EB ' { ~ <i) h{\q\) = 
TUB' 



where E c (q) = y / m 2 + q 2 . From the above expressions we have that F 2 = F 3 . 
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(B2) 
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(B4) 



(B5) 



APPENDIX C: l B ' B (\q\) AND K. B ' B (\q\) INTEGRALS 

To evaluate 1 B B (\q |) and JC B B (\q |) we use a partial wave expansion of the orbital wave functions: 

00 

1=0 

00 

*? h2 (ri,r2,r 12 )=J2fi B '(ri,r 2 )P l (n) (CI) 



1=0 
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where [i is the cosine of the angle made by r\ and r 2 and Pi(fJ-) is a Legendre polynomial of rank I. The radial 
functions /; S (Vi,r2), f t B (ri.,r 2 ), are evaluated as 

91 4- 1 f' +1 

fi B (n,r 2 ) = J i d/xP z (/x)*f fe2 (ri,r 2 ,r 12 ) 

fi B (n,r 2 ) = J i ^P(M)*f h2 (n,r 2 ,n 2 ) (C2) 

In terms of those we have 

• i B ' B (W\) 

pOO , p oo 

l B ' B (\q\) = leTr^^^UrriOOO) 2 / dr x r\ M ^ l^) / dr 2 r 2 2 ^(^fk> 2 ) 

xtffrurdtfir^n) (C3) 

with j's being spherical Bessel functions. 

• £ s ' B (k1) 

^(ki) = -^EEE E E(- 1 ) (r+ '"' +1)/ V(2L + i)(2i" + i)(2z»' + 1) 

x (H'Z" I 000) (ZV'l | 000) {I 1 LI'" | 000) : «'") 

x/ dnr^V»( -.^ klrr) / dr 2 r 2 2 |r a ) jf (n, r a ) tt L /f(n, r 2 ) (C4) 

Jo M Jo M 

with W(l"l'lL;l,l"') being a Racah coefficient and fi^ the differential operator 19 
For the actual evaluation we restrict the I, V values to I, I' = 0, • • • , 6 
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